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SPECTRAL CHARACTERIZATION OF AGING: THE REM-LIKE 

TRAP MODEL^ 

By Anton Bovier and Alessandra Faggionato 

Weierstrass Institut fiir Angewandte Analysis und Stochastik 

We review the aging phenomenon in the context of the simplest 
trap model, Bouchaud's REM-like trap model, from a spectral the- 
oretic point of view. We show that the generator of the dynamics 
of this model can be diagonalized exactly. Using this result, we de- 
rive closed expressions for correlation functions in terms of complex 
contour integrals that permit an easy investigation into their large 
time asymptotics in the thermodynamic limit. We also give a "grand 
canonical" representation of the model in terms of the Markov process 
on a Poisson point process. In this context we analyze the dynamics 
on various time scales. 

1. Introduction. The particular properties of the long term dynamics of 
many complex and/or disordered systems have been the subject of great in- 
terest in the physics and, increasingly, the mathematics community. The key 
paradigm here is the notion of aging, a notion that can be characterized in 
terms of scaling properties of suitable autocorrelation functions. Typically, 
aging can be associated to the existence of infinitely many time-scales that 
are inherently relevant to the system. In that respect, aging systems are 
distinct from metastahle systems, which are characterized by a finite num- 
ber of well separated time-scales, corresponding to the lifetimes of different 
metastable states. 

Aging systems are rather difficult to analyze, both numerically and ana- 
lytically. Most analytical results, even on the heuristic level, concern either 
the Langevin dynamics of spherical mean field spin glasses or trap models, a 
class of artificial Markov processes that in some way tries to mimic the long 
term dynamics of highly disordered systems (see, e.g., [8]). 
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One of the natural questions one is led to ask when being confronted with 
phenomena related to multiple time-scales is whether and how they can be 
related to spectral properties. This relationship has been widely investigated 
in the context of Markov processes with metastable behavior (see, e.g., [12, 
13, 14, 20, 21, 10]), and it would be rather interesting to obtain a spectral 
characterization of aging systems as well, at least in the context of Markov 
processes. To our knowledge, this problem has not been widely studied so 
far. The only papers dealing with the problem are [24], by Butaud and 
Melin, that have tackled one of the simplest trap models and on which 
we will comment below, and [17] and [23], that investigate convergence to 
equilibrium in the Random Energy Model (REM). 

The present paper is intended to make a modest step in this direction by 
analyzing the relation between spectral properties and aging rigorously in 
the REM-like trap model. While this model may seem misleadingly simple, it 
has in the past provided valuable insights into the mechanisms of aging, and 
it is our hope that the analysis presented here will provide useful guidelines 
for further investigations of more complicated models. 

The paper will be divided into two parts. In the first we analyze the REM- 
like trap model in the standard formulation of Bouchaud [9] . In the second 
part we go one step further and reformulate the model in a slightly different 
way as a Markov process on a Poisson point process. This formulation makes 
the relation to the real REM more suggestive (see [3, 4] for a full analysis), 
and allows, in a natural way, to study the dynamics of the model on different 
time scales. 

2. The REM-like trap model. Let us recall the definition of trap models 
as introduced by Bouchaud and Dean [9]. Let Q = {S,£) be a finite graph 
with vertex set, S, and edge set, £. Let E_-= ^ S} be a random field, 
called energy landscape, and let Y{t) be a continuous-time random walk on 
Q with ^-dependent transition rates, Cjj, such that Cij > iff G £, 

and 

F{Y{t + dt) = j\Y{t) =i)= dt. 

Setting '■= J^j^i'^iJ ^-iid pi^j := CijTi, the random walk, Y{t), can be 
described as follows: after reaching the site z, the walk waits an exponential 
time of mean Tj and then jumps to an adjacent site, j, with probability pij. 
In the trap model, the transition rates are assumed to satisfy the following 
properties: 

(2.1) e^»Cij=e^^Cj-i V{i,j}ef, 

(2.2) E(Ti)=oo, 

where E denotes the expectation w.r.t. the random field E_. Since in several 
physical experiments (see [26] ) the system is initially in equilibrium at a high 
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temperature, T ^ Tg and then is quickly cooled to a temperature smaller 
than Tg, and then its response to an external perturbation is measured, it 
is reasonable to consider Y{t) with uniform initial distribution. A classical 
time-time correlation function is given by 

U{t,t^) := F{Y{s) = Y{t^) Vs G [t^,t^ + t]). 

In order to observe aging, it is necessary to consider a thermodynamic limit, 
with the size of G going to infinity, and possibly a suitable time-rescaling. 
Rather recently, there have been a number of rigorous papers devoted to the 
analysis of trap models on the lattices Z [5, 18, 19] and Z"^ [6, 11]. 

In this paper we consider the simplest trap model, called the REM-like 
trap model [9], that corresponds to choosing Q to be the complete graph on 
N vertices, that is, 



N 



{l,2,...,iV}, £^:={{ij} i^jeSN}, 



and to take as energy landscape a family, E = {Ei : i € N}, of independent, 
exponentially distributed random variables, with parameter, a, with < a < 
1. Given S N, let Y]sf{t) be the continuous-time random walk on with 
transition rates Cij = e~^^/N, for i ^ j. Setting Xi = e~^', the infinitesimal 



generator of the random walk is given by 
/ (A -l)xi _xi 
N N 

X2 {N-1)X2 



(2.3) 



L 



N 



V 



N 

xn 
' N 
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X2 

' N 



N 
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The dynamics can be described as follows: after reaching the state i, the 



walk waits an exponential time of mean j^e^^ and then jumps with uni- 
form probability to another state. Although, strictly speaking, the mean 



waiting time is given by -J^e^% we call r, 



■ X, 



e * waiting time (the 



discrepancy is negligible in the thermodynamic limit A | oo). 
Note that Tj and Xi have distributions given by 



p{t) dr ■ 



■ ar 



-l-a 



dr 



(r > 1); p{x) dx 



ax 



a-l 



dx 



{0<x<l), 



respectively; in particular, E(rj) = oo. Moreover, the equilibrium measure is 
given by /ieg(0 = '^i/iJ2jLi'^j)- We are interested in the out-of-equilibrium 
dynamic with uniform initial distribution. P^r denotes the law of this random 
walk, given a realization of the random variables Ei. 

Aging in the REM-like trap model is manifest from the asymptotic be- 
havior of the time-^time correlation function 



(2.4) 



UNit,t^) ■=¥NiYNis) = YNit^) ys£[t^,t^ + t]). 
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Namely, as shown in [9], for almost all and for all > 0, 

(2.5) lim lim n7v(^t«„t«,) = ^^^^^^ ^-"(1 - nf"^ d^x. 

u,]coN]co vr Je/{i+e) 

Our main aim here is to show that the aging behavior of the system, 
derived in [9] using renewal arguments, can be obtained solely from spectral 
information about the generator L^r. The method developed below will allow 
us to get further information on YN{t) from the spectral properties of L^v- 
In particular, given a function h on (0,oo), it is possible to describe the 
asymptotic behavior of E7v(/i(2;Ar(t))) and EAr(/i(r7v(t))), where E^r denotes 
the expectation w.r.t. Ptv, and where XN{t),Tj\f{t) are defined as 

XN{t) = Xk, TN{t)=Tk \iYN{t) = k. 

These results will allow us to investigate how the walk, as time goes on, 
visits deeper and deeper traps, that is, sites with larger and larger waiting 
time Ti (see Section 2.2). 

We start by giving a complete description of the eigenvalues and eigen- 
vectors of Ljv- Let n = /ijv be the measure on Sn with ^{i) = = Ti. Note 
that Ltv is a symmetric operator on L'^{n) and, trivially, L7vl = 0, where I 
is the vector with all entries equal to 1. The following proposition is based 
on elementary linear algebra: 



Proposition 2.1. Let xi,X2, ■ ■ ■ ,xn be all distinct. Then, L^v has N 
positive simple eigenvalues = Ai <A2<---<AAr such that 

{Ai,A2,...,A,v} = {AgC:,/.(A) = 0}, 

where (/)(A) is the meromorphic function 

(2.6) ^(A):=^_l^, (AGC). 

-, X-j — A 

If the Xi are labelled such that xi < X2 < ■ ■ ■ < xi\j , then Xi < Aj+i < Xj+i, 
for i = 1,. . . ,N — 1. Moreover, for any i = 1,. . . , N , the vector i/^^*-* S M.^ , 
defined as 

^f-=^hr forj = l,...,N, 

Xj A 2 

is an eigenvector of^L^ with eigenvalue Aj. . . . form an orthogonal 
basis of L'^{fi). 

Since the Xj have an absolutely continuous distribution, we trivially have 
the following: 
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Corollary 2.2. The assertions of Proposition 2.1 hold with probability 
one for all N . 



Proof. Let A be a generic eigenvalue and let us write the corresponding 
eigenvector, ip, as ip = a(l, . . . , 1)* + u;, where J2jLi "^j = 0- Since {l^N''P)j = 
XjWj, we have to solve the system 

(2.7) = Aa + Xwj \/ j = 1, . . . , N. 

Since xi, . . . ,xn' are distinct, it must be true that a / (otherwise we get 
ip = 0). Without loss of generality, we set a = 1. Note that xj, for j = 
1, . . . ,N, since otherwise (2.7) would imply that X = = xj . Therefore, we 
get Wj = ^ . Since it must be true that J2f=i''^j = 0; we get that A is 
an eigenvalue with ip s.t. ipj = ^.^^^ being the corresponding eigenvector, iff 
(p^X) = 0. This implies that (p has at most zeros. Since (pi^O) = 0, and, for 
real A, lim^j^Xi (PW = — oo, limAfx-i 0(A) = oo, we get that 4> has exactly N 
zeros. From here the assertions of the theorem follow immediately. □ 



Proposition 2.1 has the following simple corollary: 



Corollary 2.3. With probability one, the spectral distribution '■= 
AvjL^Sxj converges weakly to the measure ax"~^ dx on [0,1]. 

Remark. The results of Proposition 2.1 are incompatible with the heuris- 
tic predictions in [24] . The discrepancy is particularly pronounced in the case 
of the eigenfunction. The reason for this is an inappropriate use of perturba- 
tion expansion in [24] . We will explain this in some detail in the Appendix. 

We will now show that Proposition 2.1 allows to derive the asymptotics of 
the autocorrelation functions easily. In fact, it contains far more information 
on the long time behavior of the systems, some of which we will bring to 
light later. 

Recall that Ptii,j), the probability to jump from i to j in an interval of 
time t, can be expressed as pt{i,j) = (e~*''"^)ij. In particular, by writing vt 
for the probability distribution of Y^v (t) and thinking of the Radon derivative 
^ as column vector, 

d^ dfi ' 



we see that 



(2 8) -a,. , 
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The above formulas are true for an arbitrary initial distribution. Taking l^q 
to be the uniform distribution, by Proposition 2.1, we get 

J N N 

^ . E where ,<^., ^ g 

Then, by Proposition 2.1 and (2.8), 

(2.9) n^(t, t.) = E E ^:^^e-«^"i)/^)^^* 
and 

(2.10) E^ihixMm = E E 

The above formulas (that may appear rather ugly at first sight) admit a 
nice complex integral representation through the following lemma: 

Lemma 2.4. Let -y be a positive oriented loop on C containing in its 
interior Xi, . . . , . Let g be a holomorphic function on a domain DcC 
with 7 C -D. Then, for any j = 1, . . . , N , 

(2-11) E = — / T7#^T dX. 

Xj — Xk 27ri (p{X){xj — A) 

Proof. Let us set X := {xi, . . . , xjy} and A := {Ai, A2, . . . , Xn}- Then, 
(/)(A) is a holomorphic function on C \ X, where i;^'(A) = X^j^i -a)^ ' and, 

in particular, (p'iXj) = ■ Moreover, the function [(j){X){xj — A)]""*^, a priori 
defined on C\ (X U A), can be analytically continued to X as a meromorphic 
function with simple poles only at the points of A. Now the conclusion follows 
from a trivial application of the residue theorem. □ 

We can obviously use Lemma 2.4 to rewrite (2.9) and (2.10) in the form 

(2.12) n„(t, t,,:) = - (av, /avj^) dx, 

(2.13) E„(K.,)) = i/,^ (av,|^/Av,^) dA. 

where Avj denotes the average over j = 1, 2, . . . , X. 

The above integral representations of nAr(t,t^) and 'EN{h{xt)) have two 
advantages. First, the appearance of averages allows to compute their lim- 
iting behavior as X | 00 easily by using the ergodicity of the random field 
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E_. Second, by means of the residue theorem, their Laplace transform can 
be easily computed in order to derive the asymptotic behavior of n7v(i,i«)) 
and K]\f{h{xt)), for A^, 1 (see Sections 2.1 and 2.2). 

A much more general derivation of the above integral representations is 
discussed in Appendix A.S. 



2.1. Aging behavior of Ilfyf{t,tyi;). 



Proposition 2.5. Let us define 



(2.14) 



1 

27d 



-tin A 



E^(e-"V(A-x)) 

A E,(l/(A-x)) 



dX, 



where is the expectation w.r.t. the measure ax" ^ dx on [0,1] and 7 is 
any positive oriented complex loop around the interval [0,1]. Then, 



(2.15) 



lim U]\f{t,tw) =U{t,t^) yt,t^ 
Nloo 



a.s. 



Proof. Recall (2.12) and fix < < 1/2. Due to analyticity, we can 
choose the integration contour, 7, to have distance 1 from the segment 
[0, 1]. For each A G 7, the random variables {xj — A)~^, j £ N, are i.i.d. and 
bounded. Therefore, for a suitable positive constant c > 0, 



(2.16) 



1 



A 



1 



A — X 



> N 



~l/2+5\ 



<e 



-cN^ 



VAG7. 



Since for each x G [0, 1] and A G 7, |^(x — A)~^| < 1, a simple chaining ar- 
gument allows to deduce from the pointwise estimate (2.16) uniform control 
in A. Using the Borel-Cantelli lemma, one can then infer that, a.s., 



(2.17) 



sup 

AG7 



A 



A 



ViVG N. 



Similar arguments show that, a.s., given M G N, there exists a constant, cm, 
such that 



sup sup 
(2.18) M-l<t<MXe'y 



-{{N-1)/N)xjt 

Xi — A 



E,. 



<cmN~^/^+' 



ViVGN. 



Note that, for each A G 7, AvjLi{xj — A)~^ is a convex combination of points 
of modulus larger or equal than 1/2, contained in an angular sector with 
angle non-larger than a suitable constant, c < vr. In particular, |Av^^(rEj — 
A)~^| > c' > 0, for all A^. From here the assertion of the proposition follows 
from Lebesgue's dominated convergence theorem. □ 
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Fig. 1. Possible integration contour. 



Given 6 > 0, we are interested in the limit of ll{0tu],tu,), as | oo. This 
will be done using the Laplace transform of Il{9ty^,tw), 



u{et^,t^)dt 

•u 



[3f?M >0]. 



The computation of this Laplace transform is trivial, if we use the integral 
expression (2.14). 

Let uj gC, with ^(uj) > 0, and fix a positive oriented loop, 7, around the 
segment [0,1], such that 7 C {z G C:?R.{z) > —?R.{lo)}; see Figure 1. Then, 
?R.{lo + X + x9)>0, for x G [0, 1] and A G 7, so that (2.14) implies 

-1 



1 

2TTi 



X{X-x){X + uj + ex)Ei 



1 



X — X 



dX]. 



Here Kx and Kx denote the expectation w.r.t. the measure qx"~^ dx on [0, 1]. 

Let us consider the change of variables z = j and write 7 for the path 7 
with inverted orientation (i.e., positive oriented w.r.t. A = cxo). Then we get 

-1 



1 
27ri 



^1 — zx){l + zuj + z9x)Kx 



1 



1 — zx 



dz]. 



Given x G [0, 1], the integrand is a meromorphic function in C \ [1, 00) that 
has only a single pole of order 1 inside 7, namely, at z = —{lo + x9)~^. By 
the residue theorem, we get 



(2.19) 



U(9,uj) 



oj + x9 + X' 



CO + xt 



UJ 



+ x9 + X 



Lemma 2.6. The r.h.s. of (2.19) is well defined and holomorphic for any 
w G C \ (—00,0]. In particular, the function II{0,lv), defined for 3ft(a;) > 0, 
can be analytically continued to the set C\ (— cx3,0]. 

Proof. As proved in [15], Chapter 3, the Laplace transform, II{9,lv), 
is holomorphic on the set of convergence points. Therefore, we only need to 
show that the r.h.s. of (2.19) is well defined and holomorphic on ^(w) 7^ 0. 
Let us assume that 9(u;) > a > 0. Then, trivially, y x,x G [0, 1], 



io + xb 



LO + x9 + X 



£B:={z£C:z= \z\e'^ with < 6* < 6*0, \z\ > c}. 
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for suitable constants, c,9q, depending on a and satisfying 6q < tt. Moreover, 
since limi^i^oo zJTMTs = 1' 



(2.20) 0<ci(a)< 



LO + xb 



uj + xO + X 



<C2(a) Va > 0, Vw : > a. 



By (2.20) and the geometry of B, we have that ^x{ J^te+x ) ^^^^ defined 
and has distance 03(0) > from the origin. Moreover, + + x| > 3^(0;). 
Therefore, the r.h.s. of (2.19) is well defined and, due to the previous esti- 
mates and Lebesgue's dominated convergence theorem, it is continuous on 
{9(tij) / 0}, thus implying continuity on C\ (— cxo,0]. 

We recall Morera's theorem: if /(w) is defined and continuous in a open 
set f2 C C and if f dio = 0, for all closed curves, 7, in U, then /(w) is 
holomorphic in fi. Therefore, using Fubini's and Morera's theorems, one 
can prove that the function ( J^q^^ ) is holomorphic on C \ (—00, 0]. The 
proof can be concluded by a second application of the same theorems. □ 

In what follows, we keep the notation, 11(0, w), for the analytic continu- 
ation of the Laplace transform. The next lemma describes the behavior of 
11(0,61;) near the origin. Using the Laplace inversion formula, we then derive 
from this result the asymptotic behavior of Ii{9,t^), as t^j | 00. 

Lemma 2.7. For any 6 > 0, set 
sin(7ra) 



A{e) =: ' ' / n-"(l - uy du. 
71" Je/(e+i) 

Moreover, define 

(2.21) ^ := {re*<^:r > 0, |(/>| <f7r}. 
Then, for a suitable positive constant, c > 0, 

(2.22) \U{e,uj)\<c\uj\~'^ VcjG^:|u;|>1, 

(2.23) \Il{9,uj)- A{9)/uj\<c\uj\~°' ^ uj £ A:\uj\ < I. 



Proof. The first estimate (2.22) follows trivially from (2.19) and (2.20). 
Let us prove (2.23) for a; € ^ and |ti;| < 1. 

In what follows, Co,ci,... denote positive constants depending only on 
6. Moreover, given z E C, we denote by Jq and the integrals over the 
paths {sz : < s < 1} and {sz : s > 1}, respectively. We extend the functions 
and defined on (0, cx)), to C \ (— cxd,0] by analytic continuation. 

Then, (2.19) implies 



LoU{e,u;) 
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X 



a-l 



X 



a-l 



[1 + Xil + e)][l + X6] 

[l + xil + 9)][l + x9r 



dy I dx 



l + x9 + y 



l/uj{l+xe) 



1 + 2/ 



dy dx. 



B := |(a;,x) G s.t. uj £ A,x = — for some s:0 < s < l|. 



Since (u;(l + x6)) ^ An {z :\z\ > cq}, we obtain 

r-l/(a;{l+a;6»)) ya-1 



(2.25) 

Jo 1 + y 

(2.26) / f— 

Let B{a) be defined as 



dy > ci, 



c^y < C2\u;{l + xe)\ 



l-a 



(2.27) 5(a) := 



oo y°'~^ 



1 + y 



dy : 



u-"'{l-u)°'-^du 



vr 



sin(7rQ) 



[note that the above second identity follows from the change of variables 
u = y(l + y)~^, while the last one is well known in the theory of the Gamma 
function]. Using (2.25) and (2.26), we obtain 



(2.28) 



Since 



< CO 



1 /•i/'^ 



dx 



B{a)Jo {1 + x{l + e)){l + x9) 



|l + x(l + ^)||l + xe|2"- 

dx 



l-a 



1/^ (l + x(l + 0))(l + xe)" 



< C4|u;| 



and, using analyticity and integrability of the singularities around z = and 
z = oo, 



x'^ ^ dx 



suj : s>0 



{l + x{l + 9)){l + x9Y 



x°' ^ dx 



{l + x{l + 9)){l + x9Y 



we get 



u;U{9,uj) 



1 



X' 



a-l 



Bia)Jo il + x{l + 9))il + x9y 



■ dx 



< C5\UJ 



l-a 
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Using the change of variables v = x ^ + 9 and u = v{\ + v) ^, we obtain 

■dx= r ^i-°(l-M)"-i(i'u, 



which imphes the assertion of the lemma. □ 

Lemma 2.7 and Proposition A.l allow us to conclude the proof of the 
aging behavior of nAr(t,t^), that is, to recover the result of Bouchaud and 
Dean [9]: 

Proposition 2.8. For almost all energy landscapes E, given 6 > 0, 

(2.29) lim lim n7v(^t^,t^) = ^^^^^ ^-"(1 - n)"-^ d^/. 

t^looNioo vr Je/{i+e) 

2.2. Visiting deeper and deeper traps. In this section we use the integral 
representation (2.13) in order to study the probability that the walk at time 
t is in a deep trap, that is, in a state with large waiting time. In Proposition 
2.9 we first prove that the probability to be in a site with waiting time smaller 
than 0(1) decays as thus implying the aging behavior of correlation 

functions described in Section 3. In the second part, we will investigate the 
random variable txN{t) and show that, for almost all E_, it has a weak limit, 
as first I oo and then i f oo. As consequence, with high probability, at time 
t, the system is in a state of waiting time 0{t), as stated in Proposition 2.10. 

Reasoning as in the proof of Proposition 2.5, we can prove, for almost 
all energy landscapes, E, that, given a function h on [0,1], that can be 
uniformly approximated by piecewise functions, 

H{t) := lim Eiv(/i(x7v(t))) 

(2.30) 

1 r e-'^J,\h{x)/{X-x))x»-Ux ^^ ^^^^ 
2TriJy A J^{l/{X-x))x'^-^dx 

where 7 is a positive oriented loop around [0, 1] . 

Since H{t) is a bounded function, the Laplace integral H{uj) := H{t)e~''^'^ dt 
is absolutely convergent when ^(uj) > 0. By the same arguments we used to 
derive (2.19), it is easy to deduce from the integral representation (2.30) 
that 

1 /o (/i(x)/(w + x))x"-idx 



(2.31) H{lv) 



to /g^(l/(u; + x))x°-idx 



In the following proposition we concentrate on the case h{x) := Ix>S- By 
(2.31), we can give precise information on the asymptotic behavior of the 
probability to be at time t in a site with waiting time smaller than 1/6: 
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Proposition 2.9. Let 

Then, for almost all energy landscapes, E, 

(2.33) lim s^-" lim Fn{xn{s) > S) = B{5)/c{a). 

sfoo A''too 

Finally, we show that, with high probability, at time t, the walk is in a trap 
of depth of order 0{t). In particular, the random variables tx]\f{t) converge 
weakly to a nonnegative random variable, as first | oo and then 1 1 oo, a.s. 
This result corresponds to the convergence of the expectation of bounded, 
continuous functions and, due to Lemma 2.11, such a convergence can be 
extended to the larger class of bounded, piecewise continuous functions, 
which is more suitable for the investigation of the phenomenon of visiting 
deeper and deeper traps: 

Proposition 2.10. Let Z be the unique random variable with range 
in (0, oo) having Laplace transform 

^, -gz-, sin(7ra) , 
E(e ) = / u [I — u) du. 

TT Je/{e+i) 

Then, for almost all energy landscape E_, given a bounded piecewise contin- 
uous function, h, on (0, oo), 

lim lim EAr(/i(ix7v(t))) 

(2.34) = lim J- / f^J.'(MxO/(A-.)K^-.i.^^ 

t^oo2mj^ X J^{l/{X-x))x'^-^dx 

= Eih{Z)). 

In particular, for almost all energy landscapes, E, 



lim lim ] 

itoo A^too 



>(^^^>n^ =P(Z<n-i) Vn>0. 



Proof of Proposition 2.9. We have to prove that lim^joo s^~°'H{s) = 
B{5)/c{a), where H is given by (2.30) with h{x) ■.= Ix>5- As in the proof 
of Lemma 2.6, we can show that the r.h.s. of (2.31) is well defined and 
holomorphic on C\ (— oo,0]. We keep the notation H for this extended 
function. Changing variables x = coy, we get 

(2.35) / ^dx = uj--' T—dy, 

Jo to + x J-i^^ + y 
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where is the oriented path {s/a;}o<s<i. Let 7^; be the path {s/a;}s>o. By 
analyticity and integrabihty of the singularities at z = 0, z = 00, we have 

,a-l 



00 y^~^ 

dy= I — — dy. 



y_ 

Let us define A := {re*^ : < r < 00, |^| < Ivr}. Then, for a suitable constant 
ci, 

,0-1 



y 



1^ 



\7u. 1 + y 



dy 



< Cl\ljJ 



l-OL 



UJ e A:\io\< 1, 



implying 
(2.36) 



1 — 1 / roo o/Q — 1 

dx = Lo^~H l^dy + O{\L0\ 
w + X \Jo l + y 



1-0- 



where A = B + 0{1/N) is understood to mean that there exists C < 00 such 
that \A- B\<C/N. Trivially, 



(2.37) 



1 rj.Ct-1 



dx 



{l + 0{\uj\)) / x'^-^dx. 



Is UJ + X 

Note that the estimate of the error terms in (2.36) and (2.37) is uniform in 
uj£A, \uj\ < 1. Then, from (2.31), (2.36) and (2.37), we get 

(2.38) |u;°iJ(u;)-S((5)| <C2|cj|i-° 



yoj eA:\oj\ < 1. 



Since, trivially, |-ff(u;)| < C3|u;| ^, for uj A with > 1, the assertion of the 
proposition follows from Proposition A.l. □ 

Proof of Proposition 2.10. As discussed before (2.30), one can show 
that, for almost all energy landscapes, E_, given a piecewise continuous func- 
tion h on (0, 00), 

<l>t{h) := lim EN{h{txN{t))) 

N'lOQ 



1 

2tH 



'tX pi 



A 



J^{h{xt)/{X-x))x''-^dx 
J^{l/{X-x))x-~^dx 



dX 



Vi > 0, 



where 7 is a positive oriented loop around [0, 1]. Note that defines a pos- 
itive linear functional on the space of continuous functions on (0,oo) that 
decay at 00 and satisfy ^>t(l) = 1. Therefore, the Riesz-Markov representa- 
tion theorem (see Theorem IV. 18 in [25]) implies that $t(/i) = Ht{h), for a 
unique Borel probability measure fit on [0,oo). In particular, there exists a 
random variable, Zt on (0, 00), such that 

lim tx]\i(t) — > Zt weakly, Vt > a.s. 
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If we take h{t) = e"*^, then $t(/i) = = U{0t,t), with IT defined as 

in (2.14). That means that Il(9t,t) is the Laplace transform of Zf. As proved 
in Section 2.1, 

l[inU{et,t) = ^^^^^ u-"(l-n)"-idu:=/(0). 

<Too TT Je/{e+i) 

We claim that f{6) is the Laplace transform of a random variable Z with 
range in [0, oo). To show this, we apply the criterion given by Theorem 1, Sec- 
tion XIII.4 in [16]. By (2.27), /(O) = 1. Moreover, f^^\e) = -2iBj[^6'-"(l + 
0)~^ , thus implying (by trivial computations) that {—l)^f^^\6) > 0. This 
completes the proof of our statement. 

Since the Laplace transform of Zt converges to the Laplace transform of Z, 
as 1 1 oo, it follows that Zt converges weakly to Z, implying (2.34), whenever 
/i is a bounded continuous functions on (0, oo). Finally, due to Theorem 5.2 
in [7], convergence still holds if /i is a bounded measurable function, whose 
set of discontinuity points has zero measure w.r.t. the distribution of Z. 
Therefore, Lemma 2.11 allows to prove (2.34) for h bounded and piecewise 
continuous. □ 

Lemma 2.11. The distribution function, F{z) ■.= ¥{Z < z) , of the posi- 
tive random variable Z is continuous. 

Proof. Trivially, F is increasing and right continuous. Therefore, it has 
a countable set of points of discontinuity. Moreover, by the Laplace inversion 
formula (see [16], XIII.4), if a; is a point of continuity, then 

F{x)= lim y iz^/{«)(a). 
^ ^ a->oo ^ n! 
'n<ax 

Given s = 0, 1, 2, . . . and 7 > 0, let 0^(7) > be such that D^a"'^ = {-lY x 
€5(7)0"'''"'^. Then the Leibniz formula implies 



(-l)"Z)^(a-"(l + a)-i) = c,(Q)c„_,(l)a-"-^(l + a)' 

< (-l)"!)^-"-^ 



-l—n+s 



Since 



/W(a) = -^^a-(l + a)-\ 
vr 



the above estimate implies 

■ / \ n— 1 

(-l)"/(")(a) = |/(")(a)| < TT (A: + a)a""-" Vn > 2. 
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In particular, given two points of continuity, < x < z, we have 



(2.39) F{z) -F{x)< — ^—^ lim sup a" 



1 



TT 



ax<n<az k=l 

One can prove that the sequence 0^=1 e~"/"(l + ^) is convergent (see [1], 
Chapter 5, Section 2.4). Denote its hmit by Cq and let 7 be Euler's constant 

7:= 



: lim 1 + - + -H 1 logn 

ntoo \ 2 6 n 



Then we can write 

n-l 



1 



n 



(2.40) 



fc=i 



n 1+ 



„o(l+l/2+-+l/(n-l)-log(n-l)) ~ 



n 



n-l 
k=l 



In particular, we can substitute in (2.40) 11^=1 e-°/'=(l + f ) with c« with 
an error term in (2.39) bounded by 

n— 1 / 



const. a "(az — ax) 



{azf 



ax 



k=l 



< c{x, z) 



jn^--(in) 



which is negligible, as a | 00. Therefore, 

F(z) — F{x) < const, lim sup CQ,a~"e'^" 

a— >oo 

(2.41) 



ax<n<az 



-l+a 



<c'(z"-x"), 

for some positive constant c'. Since (2.41) is valid almost everywhere and F 
is monotone, it follows that F is continuous. □ 

3. Other correlation functions. In this section we study the asymptotic 

behavior of different time-time correlation functions, Il^^\t,tu,), Il^^\t,tu,), 
for which deep traps play a special role. This section is mainly a preparation 
of what is to follow in the second part of the paper. 

Given 6 > 0, we define the set of sites with small waiting time as := 
{i:xi>5,i = l,..., N}. Moreover, we set 

(3.1) nS^^(t,t^) :=FN{YN{u)eDN VuG (t^,t^ + t] s.t. Yn{u) ^Yn{u-)) 



(3.2; 



FN{YN{u)£DNU{YN{t^)} 

\fU G {tyj,tyj +t] S.t. Xn{u) / Xn{u~)). 
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Given a subset A C Sjy, i £ A and s > 0, let LpN,A{h s) be defined as 

mA(^s) ■.= ^n{Yn{u) g a g [o,s]|yjv(o) = i). 

Then 

(3.3) ^ rt xe~^^^ 

+ Y.^N{YN{tu,)=i) / ds ^ ^N,Dj,{i,t-s), 
N 

(3.4) n^^(t,t^) = ^PAr(xAr(t^) =Xj)v7Ar^£,^U{j-}(j,t), 

i=i 

where the first identity can be derived by conditioning on the first jump 
performed after the waiting time and by recalhng the following realization 
of the dynamics: after arriving at the state i, the system waits an exponential 
time with parameter Xi and after that it jumps to a site in Sn with uniform 
probability. 

The following proposition is mainly a consequence of the phenomenon of 
visiting deep traps with higher and higher probability. Recall that Proposi- 
tion 2.10 implies 

(3.5) limlim¥N{xN{t)>e) = Ve>0. 

Proposition 3.1. For almost all x, 

(3.6) lim snp\U^^\t,t^)-UN{t,t^)\=0 fori = l,2. 

Proof. We consider first the case i = l. 
We claim that, for any u > and i £ Di\i, 

(3.7) ipN,DN{i,u) <exp(^-6u(^l-^-^^^^. 

In order to prove such a bound, we introduce a new random walk, Y^{t), 
whose generator, L*, defined as the r.h.s. of (2.3) with Xi replaced by 6 if 
i £ Dn- By a simple coupling argument, one gets 

where the function ^n^Djs, analogue of ipN,Df^{i,u) for the random 

walk Y^{t). At this point, it is enough to observe that ^Q^als the 

r.h.s. of (3.7). 
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Now fix e > 0. Then, due to (3.3) and (3.7), 

|nS^(t,t^)-njv(t,t^)| 

<P7v(x7v(t«,)>e) 




< PNixNitn^) >e)+e f\-Su{l~\D^\/N) 

Jo 

By the law of large numbers, 

hm re-^^(i-l^^l/^)<oo a.s. 

NIooJo 

The proposition now follows from the fact that e is arbitrary and from (3.5). 

To deal with case (ii), one proceeds in essentially the same way, decom- 
posing the path of the process at its returns to the point Xi, and summing 
over the number of these returns. One finds easily that the case when the 
process does not leave for the entire period t dominates, leading to the 
assertion of the proposition. We leave the details to the reader. □ 

4. The REM-like trap model on a Poisson point process. In this section 
we consider a slightly different formulation of the REM like trap model 
that betrays more directly its connection to the REM dynamics (see [3, 4]) 
and that offers a somewhat more natural insight in the role of time-scales 
in the analysis of aging systems. Let us consider a Poisson point process, 
V = J2i ^Ei, on R with intensity measure ae~"^ dE, where < a < 1. Note 
that such processes arise naturally as the extremal process of sequences of 
random variables. Almost surely, the support of V is an infinite set of points, 
whose maximal element is finite. Thus, we can label the points in the support 
of V in decreasing order: Ei > E2 > ■ • ■ ■ The energy landscape, E_, is defined 
as E_= {El, E2, . . .). We want to define a random process on the support 
of this point process that jumps "uniformly" from any point to any other 
point in the support. To do this, we need to introduce a cut-off. Here we fix 
an energy threshold E and set 

NE = m.ayi{i:Ei> E}. 

Note that Ne is a Poisson random variable with expectation e~'^^ . More- 
over, the probability that Ne = can be made as small as desired when E 
is chosen small enough, as we assume in what follows. 
Let Qe = {Se,£e) be the graph with 

Se.= {1,2,...,Ne}, £e:= j} -.i+j^SE). 
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Since we want to investigate the effect of time rescaling, we introduce a time 
unit, To = e^°. Then, the continuous-time random walk, YE{t), is the random 
walk on Qe having uniform initial distribution and such that, after arriving 
at site i€ £e, it waits an exponential time with mean j^^^ e^' /tq and then 
jumps with uniform probability to a different site of Se- In particular, the 
Markov generator L^; for the above defined random walk is given by L^v 
in (2.3) with := Ne and Xi := Toe'^' (since for ^ <C 0, ~ 1 when 

referring to waiting time, we disregard the coefficient j^^rj, as in Section 2). 
Note that YE{t) depends on tq, but we do not make this explicit in our 
notation. In what follows we denote by P^; the probability measure on the 
path space determined by Ye{-), and by E^; the corresponding expectation. 

Note that the physical waiting time (the absolute one) for the system at 
state i is given by Tj := e^', while in the above dynamics the waiting time is 
Tj := Tj/ro, in agreement with the choice to consider tq as our new time unit. 
In what follows we consider, when taking the thermodynamic limit E J, — oo, 
three different kinds of time rescaling: tq fixed, tq := (i.e., Eq = E) and 
To i after E [ — oo. 

As in Section 2, we are interested in the asymptotic behavior of time- 
time correlation functions. In particular, let us introduce here the correlation 
function 

IiE{t,t^) ■.= ¥e{Ye{s) = YE{t^), Vs G [t^,t^ + t]). 

We will prove that, when tq is fixed, the system exhibits fast relaxation, 
thus excluding aging behavior (see Proposition 4.2). At the other extreme, 
the scaling tq = corresponds to the implicit choice made in the standard 
Bouchaud model considered in the previous sections. In fact, with this choice 
the system can be thought of as a grand canonical version of the original 
REM-like trap model and all the results of the previous sections carry over. 
Finally, we consider the third scaling: tq | after £"1—00. In Proposition 
4.6 we show that, when performing such limits, the correlation function 
^E(t,tw) converges to f{0), where 6 = t/tw and f{0) denotes the r.h.s. of 
identity (2.29), that is, the limiting behavior of the correlation function 
II_E(i,i«;) is trivial. At this point, a simple consideration is fundamental. If 
we assume that the physical instruments in the laboratory have sensibility up 
to the time unit tq, then it is natural to disregard jumps into states whose 
physical waiting time, Tj = is much smaller than tq. Therefore, it is 
more appropriate to consider instead of HE{t^iw) the time-time correlation 

function n^^(t,t^), defined, for 5 > fixed, as 

nE^(t,t^) =¥e{xe{u) >d\/ue {tw,tw + t] :xe{u) / xe{u~)), 

where XE{t) := Xk, whenever YE{t) = k. In Section 5 we prove that U^^\t,t^) 

exhibits aging behavior: Il^^\9tw,tw) converges to f{9) after taking the (or- 
dered) limits E I —oo, tq j and tu, t oo. 
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Finally, we discuss the asymptotic spectral behavior for the above time 
rescalings. We will show that aging appears whenever the limiting spectral 
density has a singularity of order 0(x"~^) at 0. 

Let us recall some properties of the Ppp with intensity measure 

aTQ°'x°'~^ dx on (0, oo), which will be frequently used below. Given M > 0, 
the truncated Ppp, J2xi<M ^Xi, can be realized as follows: Let um be a Pois- 
son variable with expectation (^)" = aTQ°'x°~^ dx and let Xi, z G N, 
be i.i.d. random variables on [0, M] with probability distribution p{X) dX = 
aM-°X°-idX. Then 

riM 

(4.1) E^-»~E^^- 

Xi<M i=l 

in the sense that the point processes above have the same distribution. In 
particular, taking M = Toe~^ , we get 

(4.2) E^-»~E^^- 

i<NE i=l 

where is a Poisson variable with expectation e""^^, and G N, are 
i.i.d. random variables, independent of n^, distributed on [0,roe~^] with 
probability distribution p{X) dX = e"^aro""X"-i dX. 

Notation. It is convenient to introduce the random walks XE{t), TE{t), 
defined as 

XE{t):=Xk, TE(t):=Tk iiYE{t) = k. 

We denote by je the positive oriented loop having support 

supp(7£;) = {x±i:x£ [-l,Toe~^ + 1]} 

U {-1 + bi : \b\ < 1} U {tqc-^ + 1 + 6i : |6| < 1}. 

Moreover, we call 7oo the infinite open path, oriented from oo + i to oo — i, 
having support 

supp(7oo) = {x±i:x> -1} U {-1 + bi : \b\ < 1}. 

Finally, for given E, = x[^^ < Ag^-* < • • • < A^^^ are the Ne distinct 
eigenvalues of the infinitesimal generator L^; (see Proposition 2.1). 
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4.1. To fixed. Let us first observe that J^i^i'^i < °o, for almost all E. In 
fact, since the Ppp J2i^Ti has intensity measure aTf)T~^^^"^ dr on (0, oo), 

/oo 

e( y Ti]= /\ro"r-°dr<oo. 

Whenever J2'^iTi < c>0; it is easy to derive the asymptotic spectral be- 
havior of the system from Proposition 2.1 and to show its fast relaxation, 
thus implying the absence of aging: 

Proposition 4.1. For almost all E, 

oo 

(4.3) hni^ y (5^{E) = y Sx^ vaguely in A^([0, oo)), 

where A^([0, oo)) denotes the space of locally bounded measure on [0,oo) and 
{0 = Ai<A2<A3<---} = |AGC:f;-^ = o|. 



k: 



Proof. In what follows we assume that X^i^i Ti < which is true a.s. 
Then the function (f>oo{X) := X^fcLi x^-\ i^ ^^^^ defined on C \ {xi :i>l} and 
has nonnegative zeros = Ai < A2 < • • • , such that < A,/ < Xi for any 
i> 1. At this point it is enough to show that 

lim Ap^ = Ai V« = l,2,.... 

Ei-00 

The assertion is trivial for i = l. Suppose that A'^^; > i > 1 and set iPe{X) := 

]k=i Due to Proposition 1 ^^^^ 

interval (xi_i,Xj). In particular. 



J2k=i Xk~\ • t° Proposition 2.1, A,- is the unique zero of iPe{X) in the 



i 

Since ipEiX) > for ah A G (xj_i,Xj), we get 

\xf^ -Xi\<{xi-x,.i)^\i^E{Xi)\ 

and, therefore, the assertion follows by observing that the identity (pooiXi) = 
implies 

00 -j^ 

\i^E{Xi)\< -^—^^ as^i-00. 

k=NE+l n 
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Proposition 4.2. For almost all E, 



(4.4) lim lim ¥E{xE{t) = xj) = — ^ 

tloo El-oo l^k=l^k 

thus implying 



(4.5) 
(4.6) 



lim lim TlE{Otw,tw) = ^ 

iu)too El—oo 



lim lim n^(t,ti 



Xjt 



Vj = l,2,... 

V6'>0, 
Vt >0. 



Proof. In what follows we assmne that E_ satisfies Yli'^i < oo. Setting 
h{x) = lx=xj in (2.13), we get the integral representation 

-1 

d\. 



(4.7) FE{xE{t) = x,) = ^.j 



-\t 



N 



^ 1 



IE K^j ~ \ i^^i — ^ 

Applying the residue theorem [see the arguments used in order to derive 
(2.19)], it is easy to compute the Laplace transform, Fe{oj) = /g°° WE{xE{t) = 
Xj)e~'^^dt for 3f?(cj) > 0: 

(4.8) Fe{u:)=Uu + x^)Y.-—-] ■ 



Using that, ioi uj = a + ih and N is any positive integer 

1 

1 ' ' 



N 

^-.UJ + Xj 

J=l J 



> 



\U! + Xl\ 

b 



if a > 0, 
if a < 0, 



(a + xi)2 + 62' 

we obtain that, almost surely, there exists c> 0, such that 



(4.9) \FE{io)\<o 



1 



yE, Vcj G ^ := {re^^ : < r < oo, |e| < |^}. 



\L0\ 



Let us now introduce the path, 7, consisting of the parabolic arcs {— i ± 
it"^ :t> 1} and the circular arc of radius around the origin connecting 
(in anti-clockwise way) —1 — i to —1 + i. The orientation of 7 is such that 
— 1 + i comes before —1 + i. Then, by means of (4.9), the Laplace inversion 
formula and Lebesgue's dominated convergence theorem, we get 

^lim ¥E{xE{t)=Xj) — — 



El-QO 



2m 



e"^F{Lo)du;, 



F{uj) := u;{oJ + x^) ^ 



k=l 



UJ + Xk 
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Note that F{u>) is the hmit of Fe{uj), as E I oo; in particular, it satisfies 
(4.9). Moreover, F{ui) is the Laplace transform of limEi-oo^EixEit) = Xj) 
and 



ujFiu) 



< cIcjI V IwI < 1 : u; S 



At this point (4.4) follows from Proposition A.l. Moreover, from (4.4) and 
the identity 

Ne 

UE(,t,t^) =Y.FE{xE{t^) =xj)e~^''^-'y''^^^\ 
i=i 

one infers (4.5) and (4.6). □ 

4.2. To = e^. Note that, choosing To = e-^, the random variables Xi, . . 
introduced in (4.2) are i.i.d. with distribution given by p{X) dX = aX"'^ dX 
on [0, 1]. Therefore, due to (4.2), we can think of 1^(0 as the grand canon- 
ical version of Bouchaud's REM-like trap model. In particular, it exhibits 
the same asymptotic spectral density and the same aging behavior: 

Proposition 4.3. For almost all E_, 

Ne 



\\m. (5^(E) = ax" ^ dx weakly in j^A {[0,1]). 



Proof. Approximating continuous functions on [0, 1] by step functions 
having rational values and jumps at rational points, it is enough to prove 
that, given < a < & < 1, 

hm ^ I {j : 1 < j < AT^ , Af ) g [a, 6] } | = 6" - a" a.s. 

EJ.-00 1\e 

We set 

AE:=\{j:l<j<NE,Xje[a,b]}\ 

= |{j:j>l,e-^^G[e-^a,e-^6]}|. 
Then, due to Proposition 2.1, we only need to prove that 

lim — — = b°' — a.s. 

Ei-oo Ne 

To do so, observe that Ne and Ae are Poisson variables with parameter 
e""^ and e-"^(6"-a"). For n G N, we set E{n) = Inn, that is, e-°^(") = 
n^. The Chebyshev inequality and Borel-Cantelli lemma then imply that 
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By monotonicity, one can extend the first limit to lim£;|_oo J^^e = 1, a.s. In 
order to extend the second limit to general we observe that, whenever 
E{n + l) <E<E{n), 

<|{i:e-^^- G [ae-^("),ae-^("+i)]U[6e-^("),6e-^("+i)]}|. 

Since the r.h.s. is a Poisson variable with expectation of order 0{n), the 
Chebyshev inequality and the Borel-Cantelli lemma imply 

lim sup — ~r~pv~r^ ~ 

n1ooE(n+l)<E<E{n) e 

which allows, to prove that lim£;|_oo ^-aE =b" — a°, a.s. □ 

Proposition 4.4. For almost all E, and 7 a positive oriented loop 
around [0, 1] , 

lim UE{t,tu]) 
Ei-00 

(4.10) ^ /■ e-*-^ /o (e'^'VlA - x))x''-^ dx 



2TTiJ^ A /oi(l/(A-x))2;"-i(ix 
In particular, for almost all E, given 9 > 0, 

(4.11) lim lim UE{et^,t^) = ^^^^^^ r tx-"(l-n)"-^du. 

t^Too -Ei-00 vr Je/ii+e) 

Proof. Our starting point is (4.2) and the following inequality, which 
holds for any bounded function, /, with E(/(Xj)) =0: 

P(|Av^=i/(Xj)| >5) <2exp y5>0,k = l,2,.... 



In particular, conditioning on n'^ [see (4.2)], we get 

(4.12) P(|Av"f J(Xj)| > 6) < 2exp{-e-"^(l - e"'^'/^"^"- )}. 

(4.10) can now be derived from (4.12), the Borel-Cantelli lemma, and the 
integral representation 



1 /• e-*»^Av:^^i 6-^^7(3;. -A) 

(4.13) UE{t,t^) = — / — N,,^ ' d\ 

27njy A Av^ J^l/(xj - A) 



where 7 is a positive oriented closed path around [0,1] [see (2.12)]. Note 
that the r.h.s. of (4.10) corresponds to the function n(i,i^) introduced in 
Proposition 2.5. Therefore, the assertion of Proposition 4.4 follows from 
Propositions 2.5 and 2.8. □ 
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4.3. To i after E [ —oo. In this scaling regime, we show that the vague 
hmit of the suitably rescaled spectral density is given by the measure ax"~^ dx 
on [0, oo) and we recover the aging property of the correlation function. 
Moreover, due to the fact that we are effectively already at "infinite times" 
on the microscopic scale, we get a pure aging function even before taking t 
and tw to infinity: 

Proposition 4.5. For almost all E, 

Ne 

lim lim Tq 6 (e) = ax'^'^ dx vaguely in ^A{[0,oo)). 

TolOEl-oo ^ 

Proposition 4.6. For almost all energy landscapes, E_, given positive 

t, ty^ , 

4.14 hm Hij t,t^ =— / — ^'oo 1// T^dX 

and 

lim lim IlEit.ty.) 

TolOEl-oo 



^'•''^ u-{l-ur-Uu u,heree=L 

vr Je/{i+e) t^ 

Remark. The integral in (4.14) exists due to Lemma 4.8. 



Due to Proposition 2.1, Proposition 4.5 follows, if one is able to prove 
that TQj2f=i^Xj converges vaguely to ax°'~^ dx on [0,oo) when taking the 
(ordered) limits E | — oo, tq | 0. This is the content of Lemma 4.7 below 
(which is analogous to Lemma 4.16 in [4]). Finally, the proof of Proposition 
4.6 is based on (and given after) the technical Lemmas 4.7 and 4.8. 

Lemma 4.7. Let M > and let f be a bounded, continuous function 
on [0,M]. Then, there exists 5 > 0, such that, for almost all energy land- 
scapes, E, 

rM 

ro E /(^i) - / /(^)«^""' dx 



(4.16) 

Xi<M 

where c > is a positive constant. 



<CT^ Vro>0, 



Proof. Let Xi,X2,-.- and um be as in (4.1). Due to (4.1) and since 
Var(nM)=E(nM) = (M/To)", 

\{j:xj<M}\ 



(M/to)° 
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In particular, given 7, s > such that 2s — 7a < — 1, using the Borel-Cantehi 
lemma, we obtain that, for almost all energy landscapes, there exists c > 0, 
such that 

\{j:Xj<M]\ ^ 

Due to this estimate, we get that 



< ck ^ V A; = 1, 2, . . . , where To := A; 



(4.17) 



To E fixi)-M^Ay,^<Mfix^, 



Xi<M 



< ck^ 



V/c G 1, 2, . . . , where To := A; 



where Avxi<M denotes the average over the set {xi < M}. As done for (4.12), 
if 0<p< 1, 

rM 



Av,,<M/(xi) - / /(x)ax°-^ dx 

Jo 

S2e.p{-(M)'"a_e--,} 



> 



In particular, by the Borel-Cantelli lemma, for almost all E, 

rM 

' < ck- 



(4.18) 



Av,,<m/(xO - M— / f{x)ax''-' dx 
Jo 



Vfc = 1, 2, . . . , where To := /c 

if s is chosen small enough. Now (4.17) and (4.18) imply the assertion of the 
lemma, if tq = k~"' , for some A: = 1, 2, ... . The general case tq > follows 
easily from the uniform continuity of /. □ 



Lemma 4.8. For almost all energy landscapes, E_, there are positive con- 
stants, rQ,ci,C2, such that the following holds: If tq <Tq , N > \{j -.xj < 1}| 
and A G 700 (or \ = a + ib, with |6| < 1 and a>Xj + 1, for all j < N), then 

N 



Moreover, if tq < Tq and M > 1, then 



roE 



1 



>ci\X\-\ 



(4.20) ro"E 



1 



1 \xj - M 



<C2|Ar"Mn(l + |A|) */AG7o 



(4.21)ro" E I — ^T<C2M"-MnM i/ A G 700 or 3fi(A) = M + 1, 
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where Y.x,>M means Ej>i:x,>M- 

Proof. It is convenient to introduce the nonrescaled Ppp J2i ^y^i where 
yi:=e~^\ with intensity measure ay°'~^ dy on (0,oo). We set xj = TQyj. 
Moreover, we fix /3 > 2 and 0<7</3/2 — 1 and we define 

iVn:=|{j:n^/"<y,<(n+l)^/"}|, 

for n positive integer. Then the Borel-Cantelh lemma imphes that, for al- 
most all E_, 



(4.22) 



iVr, 



1 



<cn^ Vn = 1,2, 



Moreover, again using the Borel-Cantelli lemma and a simple argument 
based on monotonicity, one can prove that there exists (5 > 0, such that, for 
almost all E, 



(4.23) 



|{j:2/j <n}| 



1 



< KU 



V'u> 1, 



where k is a positive constant. We leave this proof as an exercise. 
In what follows we write \ = a + ih. Then 



(4.24) 



N 



1 



^ N 

E 



N 



+ 6^ 



In order to prove (4.19), we assume (4.22) and (4.23) to be valid and 
we let < To < Tq < 1, where Tq is such that 1 > k{tqY . This implies that 
{xjixj < 1} /0. By (4.24), ifN>\{j:xj < 1}\ and AG 700, 



TV ^ 

-, Xn A 

3=1 J 



> < 



E 



{xj — a)2 + 6^ 



>c|A|-Vo°|{x,:x,<l}|, if |5|>i. 



^oE (,^.:\;V,. >cro1{x.:x,<l}|, 



if |6|<i. 



At this point, (4.19), for A € 700, follows from (4.23). The case \ = a + ih, 
with |6| < 1 and a > + 1, for all j < N, can be treated similarly. 

It is easy to derive (4.20) and (4.21) from the estimates (4.25)-(4.31) be- 
low that hold for almost all E: 
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if a < 100, 1 < M and X = a + ib£ 7oo, then 



(4.25) 
(4.26) 



1 



fr{ \xj - A| 



,>A/ I 



Xj — A| 



< cM' 



a-l. 



if To is small enough and a > 100, then 



(4.27) 



ro E 



(4.28) To" 

a/2<Xj <a—l 

(4.29) To" 5: 

a— l<a::j<a+l ' 

(4-30) E uT^ 



(4.31) 



-o"E 



1 






a\ 




1 








a\ 




1 








M 




1 




\Xj 




a\ 




1 




\Xj 




a\ 



< ca"~^ In a, 

< ca"~^ if A = a + i6 G 7o 



< ca" ^Ina, 



< cT 



a-l 



if T > 2a. 



Let A e 7oo with a < 100. Then, due to (4.23), 



-0 E 



j ■ Vj < 



To 



<c', 



while, due to (4.22), 

(4-32)ro"E^^<-o"E:^<-'-o"^ E n^"^"^/" < c", 



x,>l 



-.1 -"/|3 I 



which proves (4.25). The proof of (4.26) follows from the same arguments. 
(4.27) is a simple consequence of (4.23). The l.h.s. of (4.29) can be bounded 
hy TqIIJ : a — 1 < tqUj < a + 1}| and (4.22) allows to conclude the proof of 
(4.29). 

The proof of (4.28), (4.30) and (4.31) can be easily derived from the 
following estimate. Let 1<A<B with S<a-lorA>a + l, then (4.22) 
implies 



^0 E r^<c^o T r 



u \a — Tox 



/3/« 



■ dx 
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where n_ = L(^/to)"/^J - 1, n+ = L(^/^o)"/^J + 1, n = n„ - 1, = n+ + 1 
(we assume tq small enough in order to exclude the singular point in the 
above intervals of sum and integral). □ 

Proof of Proposition 4.6. In order to avoid confusion, we underhne 
here the dependence on tq = e'^" by writing IlE,Eo{tj'l^w) instead of n£;(t,t^). 
Our starting point is given by the integral representation (2.12): 

(4.33) UEMt,t^) = l^ , a^k,, dX. 

Let us choose E satisfying Lemma 4.8. Then, due to the exponential decaying 

T* 



factor e and Lemma 4.8, if tq < Tq , and E is small enough, such that 



tqc"^ > 1, the integration path 7^; in (4.33) can be replaced by 700. At this 
point, (4.14) follows from Lemma 4.8 and Lebesgue's dominated convergence 
theorem. 

To prove (4.15), given a positive integer M, we set 

,,_1 f e-*-^ro"E.,<Me-"^V(^i-A)^^ 
9MMt, t^) - ^ A ro"E.,<Ml/(x.-A) 

where Tm is the positive oriented path whose support is 

supp(rj\,/) = {A G C : |A — x| = 1 for some x G [0,M]}. 
Then, applying Lemma 4.8, whenever tq < Tq , 



lim IlE,Eo{t,tw) - gM,Eo{t,tw) 
El— 00 



< cM"-^ In M VMeNu 



Let us assume that E_ satisfies (4.16), for all M £ N4. and for f{x) = ^rrx: 
or f{x) = jta' foi' A in a countable dense set of Fm and for all rational 
positive t. Then, by a chaining argument, we get 

lim gM,Eo{'t,tw) = gM{t,tw) yt,tw>0, 

Eol-oo 

where 

.X . N 1 /■ e-*-^/(f^e-^V(A-x)x°-idx „ 

4.34 g^,{t,t^)= rM.u^ , „ , , dX, 

2mJTM A l/(A-x)a;°-i|ix 



which implies 



lim sup 

-Eoi-00 



lim TVE,Eo{t,tw) - guit^tw) 
E[-co 
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At this point, we observe (see the proof of Lemma 4.8) that there exist 
c, c' > such that for all M G N+: 



r-M 




-1 


Jo 


A- 


X 






-1 


10 


|A- 


X 


POO 


x"- 


-1 


IM 


|A- 


X 



dx 



>c\\\-^ VAerA/U7, 



OO ) 



(4.35) / rdx<c! VAsTmUtoo, 



OO- 



From the above estimates, we infer 
(4.36) |(7M(t,t»)-5(t,i«-)l <cM°-ilnM, 



where 



e 



(4.37) ^ j^^ — _ ^A. 

Using the analytic properties of the integrand in the r.h.s. of (4.37), one 
can show that g{t,tuo) = git/tw, 1). In order to compute g{9, 1), we observe 
that for a suitable positive constant c, \g{9s,s) — gM{ds,s)\ < cM"~^lnM, 
for any s > 1 [in fact, the constant c in (4.36) can be chosen uniformly if 
tw > 1]- By the results of Section 2.1 [cf. (4.34) with n(t,i 

w) Proposition 

2.5], we get 

lim gM^ds, s) = r.h.s. of (4.15), 

OO 

thus concluding the proof. □ 

5. Other correlation function when tq J, after E I —oc. As stated in 
Proposition 4.6, the standard time-time correlation function IlE{t,tw) has 
trivial behavior after taking the limits ii^ | — oo, tq i 0. For physical reasons, 
it is more natural to disregard jumps into states with physical waiting time 
Tj = e^"- much smaller than tq , since we assume that the physical instruments 
in the laboratory have sensibility up to the time unit tq. Therefore, let us 
fix (5 > and consider here the more natural time-time correlation function 

{t, t^) = ¥E{xEiu) >Syu£ {t^, t^ + t]: xe{u) + xe[u~)). 

The main result of this section is the following: 

Proposition 5.1. For almost all E, 
(5.1) lim suplimsuplimsup |n^^ (t, t^„) — n£;(t, t^)| = 0. 

tw^o^ t>0 Eoi-oo Ei-oo 
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In particular, for almost all E_, 
(5.2) 



lim lim lim Yi^^ {9tw,t^ 
tw]oo Eal—oo Ei~ao 



sin(7ra) 



Je/a+e) 



TT Je/{i+e) 

Note that the correlation function Il^^\t,t^) is the analog of nj^^(t,t^) 
of Proposition 3.1. As for the proof of Proposition 3.1, a useful observation 
is that, given 5 > 0, 

(5.3) lim lim lim FE{xE{t) > 6) =0 a.s. 

We can prove a stronger result concerning the phenomenon that with high 
probability the system visits deeper and deeper traps. In fact, note that, 
by (2.13), 

> ^) = 2^ L ^^^NTTV, —dX. 



7b 



^ E?Jil/(x,-A) 



Then, reasoning as in the proof of Proposition 4.6 and using the results of 
Section 2.2, one can easily show the analog of Proposition 2.9: 

Proposition 5.2. For almost all E, 



(5.4) lim lim lim t'-"'¥E(xE(t) > 6) - , 

nooT-oiOEi-oo c(a) 

where 

foo a-2 ^ poo 

B{S):= rooZ-ifnl\^ ' c(a):=/ y^-^e-^dy. 
Jo x° ^/{l + x)dx Jo 

Proof of Proposition 5.1. Trivially, (5.2) follows from (5.1) and 
Proposition 4.6. Since Eq | — oo after E [ — oo, we assume that E < Eq and 
define 

DE,Eo-={i-E<Ei<Eo}. 

This set corresponds to the small traps into which we allow the particle to 
jump. By (5.3) and the same arguments used in the proof of Proposition 3.1 

for ? = 1 (with exclusion of the last step since here lim^|_oo ^^pj'^^^ = 1 a.s.), 
it is easy to derive the assertion of the proposition from Lemma 5.3. □ 
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Lemma 5.3. For almost all E_, there exist positive constants, p,c, in- 
dependent of E,Eq, satisfying the following property. Whenever \{i:Ei > 
^o}|>0, 

(5.5) limsup— V (pE,Eoiht)<ce-P^ Vt>0, 



■ E, 







where De,Eo '■= {i'-E < Ei < Eq}, for E < Eq, and the function (pE,Eo is 
defined as 

(5.6) ^E,Eoii,t) ■.= ¥EiYE{u)£DE,Eo yue[0,t]\YEiO) = i). 

Proof. Let us assume that E < Eq, \{i:Ei > Eq}\ > and, without loss 
of generahty, 6 = 1. 

We fix £ > such that e~°^ < ^ and define 

Wi,E:={i:E<E^<E + e}, iVi,^ := \Wi,eI 
W2,E ■.= {i:E + e<E,<Eo}, N2,e := \W2,e\. 

Note that De,Eo = Wi^e U W2,e and Ne, Ni^e, N2,E are Poisson variables 
with parameters e~°^^ , e~°^(l — e~"^) and er°^^~°^^ — e""^", respectively. 
In particular, for almost all E_, 

lim ^ = 1, pi := lim ^ ' = 1 — e , 

(5.7) 

r N2,E -ai 1 

P2 := lim — =e ""^ < -. 

Ei-oo Ne 2 

We observe that h := Ne — Ni^e — I^2,E is a positive integer, independent of 
E and Xi > e.^o-E-i ^ if i g Wi^e, while > 1, if i S W2,e- Let us introduce 
a new random walk, Y^{t), on Se, whose infinitesimal generator, L|;, is 
defined as L^;, with Xi replaced by x* defined as 

Xi, \ii^DE,Eo, 
Ae := e^o-^-^ if i E Wi,E, 
1, if i E W2,E- 

We denote by the probability on path space associated to Y^{t) with 
uniform initial distribution and we set 

(5.8) ^*E,Eoii,t) :=P^,(Y^(n) E De,Eo Vn E [0,t]|yi(0) = i). 

By a simple coupling argument, one gets (fE^Eoiht) ^ ¥'£;,Eo(^' partic- 
ular, 

(5.9) ^ iei^B.Eo ieDe.B, 

yieDE,Eo, Vt>0. 
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At this point, it remains to estimate In order to simplify notation, 
we write D, N, Ni, N2, A, dropping the index E. Moreover, we consider the 
following realization of the dynamics of Y^: after arriving at a site i, the walk 
waits an exponential time of parameter x* and then it jumps to a point of 
Se with uniform probability. In particular, jumps can be degenerate, that 
is, initial and final sites can coincide. 

We claim that 



where 



(5.10) 







*2 := 2^ 2^ 



(5.11) 



A:i=0 k2=l 



N J \N 



fci! (K2-I)! 



(5.12) $3 ■•=^e.,{-At(l - f ) } + f exp{-«(l - f 

The above identities can be derived from the probabilistic interpretation of 

ki, k2 as 

ki = I {jumps performed before time t having starting point in Wi}\ 
and from the following simple identities: 
P(ri +T2 + --- + Tn£[z,z + dz)) 

^n— 1 

= e K -. -T-j- dz, 
[n — Ij! 

P(ri + T2 + ■■■ + Tn < z and Ti + T2 + ■■■ + Tn + T^+i > z) 
n\ 

where z > and Ti,T2,. . . are independent exponential variables with pa- 
rameter K. 
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Finally, we only need to prove that, for suitable positive constant c,p> 0, 
(5.13) limsup^i <ce~P* Vi = 1,2,3. 

-Ei-oo 

We give the proof in the case i = 1; the case i = 2 is completely similar, 
while the case i = 3 follows directly from (5.7). 
We fix 7 : a < 7 < 1, set ko := A'^ and write 

where ^y^° is the contribution to of addenda in the r.h.s. of (5.10) with 
1 < ^1 < ^0 and k2 > 0. 
If ki > ko, then 



Nj \ N J \ N-N2J -\ N J \ N 
thus implying that 

< (^,A..) (1 - 1)^' < (1 - 1)^' i ^El-^. 

where ^i(^^^,iV2) is defined as in the r.h.s. of (5.10) with A'^i replaced by 
N — N2. Note that it does not exceed 1 since it corresponds to the probability 
of a certain event. 

Let us now consider the term <I>^^". To this aim, since {^^^^^) < 2'^'^'^^', 



^<ko^ ^ f2Am\''^f2tN2\''' l{0,t ) 



ki=lk2=0 

where 



)!fc2!' 



I UJ2 



Fix < m < 1 with 2p2 + 2mpi < 1 (recall that 2p2 < 1). Then, trivially, 
,5.14) ,(„-), 

From such a bound, one gets immediately 

^ ^ f 2ANi \''U2tN2\''^ 1 / nU 

tito^ M iV J {k^-iy.k2l [ ' A 

(5.15) 

' ' m Ni N2 



<cexp|-t^l---2-n.-2- 

The last expression, when E [ —00, converges to cexp(— i(l — 2pim — 2p2)), 
in agreement with (5.13). 
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In order to estimate the integral I{^,t) = e~* J^t/A e~^^~^'^'^u'^^~^ du, we 
observe that 

e-'Vdu = {-iy 



dz^ z 

< hs+- H iw + -] ys,w,z>0, 



Z \ Z J z \ z 

thus implying the bound 

The contribution of ce~^A-^^{mt + l)''^~^ to ^>f''" can be treated by means 
of estimates similar to the ones leading to (5.15). 
In order to conclude, we only need to show that 

^ ^ f2Am\''~'f2m\''^__ I „ 

(5.16) .^i.^A ^ ^ V iv ^ {h-m,i^ 

as £' I — oo. 

To this aim, observe that 

r.h.s. of (5.16) < e-^*+2t^^/^ |^ {^^^ 

< c{t)e-^^A^{4tA)^'' =c{t) exp{- At + -fin A + A^ ln{U A)}. 
Since < 7 < 1, we get (5.16), thus concluding the proof. □ 

APPENDIX 

A.l. Laplace transform. 

Proposition A.l. Let G{t) be a bounded measurable function on (0, 00) 
and let us consider the Laplace transform 



G{uj)= / G{t)e~*'^dt 
Jo 



well defined if 3?(w) > 0. Let us define 

(A.l) ^:={re'^:0<r<cx),|(9| <|7r}. 

Suppose that G can be analytically continued to C \ (— oo,0] and that there 
are positive constants 7, /3, a, c and B gM. such that 

(A.2) \G{io)\<c\uj\-^' VcuG^,|cj| > 1, 

(A.3) |w^G'(cj)-S| <c|cj|" VcugA|cj|<1. 
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Then, 

lim s^-f^Gis) = — — where c((3) := / /"^e"^ dy. 

Proof. If we set H{s) := ^rpys'^"^ with s > 0, then the Laplace trans- 
form H{oj) is well defined for ^{uj) > 0, H{uj) = Blo~^ and, trivially, H{uj) 
can be analytically continued to C \ (— oo,0]. 

By the inverse formula for Laplace transform (see Chapter 4, Section 4 
in [15]), we have 

(A.4) G(s)=lim— / e''^G{iu)dLO Vs>0,x>0, 

K-^oo 2m Jx-iK 

where w runs over the vertical path connecting x — iK and x + iK. The 
above formula remains true if substituting G with H. Therefore, 

R i-x+iK 

s^~^Gis) - ^ = lim / e'^iG{uj) - H{u;)) duj 

Vs>0,x>0. 

Let p := min(7,/3)/2. Fix a positive number x and, given K and s, define 
the following paths (see Figure 2). 

'jK is the vertical path from x — iK to x + iK. 7i^_|_ is the segment from 
—s~^ + s~^i to —1 + i. 72,+ is an arc from — 1 + i to —K^ + iK given by 
the parametrization z{t) = —t + it^/P with t G [1, KP\. 73,+ is the horizontal 
segment from —K^ + iK to x + iK. For i = 1, 2, 3, we define the path 7j,_ 
by considering the reflection of 7j,+ w.r.t. the real axis and inverting the 
orientation. Let 70 be the positive-oriented circular arc of radius from 
— s~^ — s~^i to —s~^ + s~^i crossing the axis of positive real numbers. 

Note that the above paths depend on s and/or K. 

Because of analyticity, the integral over of e^^G{uj) is equal to the 
sum of the integrals over 73,-, 72,-, 7i,-) 70; 7i,+ ; 72,+ , 73,+ - The same is 
valid with G replaced with H. 

By (A. 2), we have that 

(A.6) / \e'''G{uj)\\duj\<ce'''KP-^ [{) asK^oo 

•^73, ± 

and, for a suitable rational function /, 
^1-/3 j \f,s^G{uj)\\duj\ 

"'72,± 

/oo , 
^^s{-t+it ^^G^-t + it^/p){-l + p'H^^-P^/Pi)\ dt 

/oo 
e-'^f{t)dt<c's^-'^e-'^^ 10 ass Too. 
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-Rf' + iK 



X + iK 




Fig. 2. The integration path in the Laplace inversion formula. 



Similarly, it can be proved that the corresponding integrals with G substi- 
tuted with go to by taking the limits K ^ oo, s j oo. 

Let us now estimate s^~^ ^^-i e~^^t°'~^ dt by dividing the path of integra- 
tion in two paths. Choosing < 5 < 1, 

-l+S 

,1-13 r ^-st^ 



.'-a _|_ „'~a+5{l+a-f3) 



< CS-" + CS 



a-/3 



-stj.a—f3 



e-^ dt < s 



9{s) 



for a suitable rational function g{s). In particular, choosing 5 small enough, 
the above upper bounds imply 

lims^-^ e-'H''-f^dt = 0. 

stoo 

This result, together with (A. 3), implies 

^ / e''^{G{uj)-H{uj))duj = 0. 

ni,± 



lims^"^^ . 

sjoo ZTTl 
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Trivially, by (A.3), 



,1-/3 



70 



<s°|0 asst oo- 



The proposition now follows from the estimates above and (A. 5). □ 



A. 2. Perturbation theory. In this appendix we comment on a paper by 
Melin and Butaud [24] where the eigenvalues and eigenfunctions of the gen- 
erator of our model were computed using perturbation theory. As pointed 
out earlier, these results are at variance with our exact results, and it may 
be worthwhile to point out the flaw in their arguments. Melin and Butaud 
write the generator L defined in (2.3) as L = T + j^T^^\ where 



T:-- 



( x\ 


Vo 

/ 





X2 



\ 




Xl 
X2 



\-XN 



-Xl 
-X2 

-XN 



-Xl ^ 

-X2 

-XN ) 



The factor in front of the second term encourages them to consider this 
term as a small perturbation. Both T and T^^^ are symmetric operators on 



L ifi), where //(i) 



We denote (•,•) the scalar product in L (fj,) and 



assume xi, . . . ,xn to be distinct positive numbers. 

Given an operator A : — > we write || A|| for its operator norm. 

Because of symmetry, ||r|| and are given by the maximum of |A|, 

with A eigenvalue. Trivially, T has eigenvalues xi, . . . ,xn and T{ei) = Xiei, 
where ei,...,eM is the canonical basis of M^, while T^^^ has eigenvalues 
0, -{Xi + X2 H h xn). 

Given z E C, we can define the holomorphic function T{z) = T + zT^^\ 
A natural condition in order to apply perturbation theory to T{z) (see [22], 
Chapter II) is 



(A.9) 
where 



z < 



d 
2ao' 



d = inf \xi — X 



On := min ||T' ' — a\ 

am 



Xl+X2-\ \-XN 
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In this case, we can conclude that T(z) = T + zT^^^ has eigenvalues 
Ai(z), . . . , Xn{z) with Afc(z) = j:n=o 4"^^''' where, for n>2, |A^"^| < a^{2/dY'-\ 



and 



^(1) ^ (r(^)efc,efc) 
(e/c,efc) 



{ek,ek){ej,ej) ' 



Similar series exist for the perturbed eigenvectors. 

However, the crucial condition (A. 9) is hardly satisfied when z = j^, since 
it reads 

(A. 10) Avf^iXj < inf \xi - Xj\, 

while a.s. the l.h.s. of (A. 10) has nonzero limit and the r.h.s. converges to 
at least like 1/A^. 

The fact that the conditions for the application of perturbation theory 
are violated explains why its predictions are incorrect. This discrepancy 
happens not to be too obvious as far as the eigenvalues are concerned (which 
are caught between the diagonal elements of the generator and thus are 
somewhat similar to them, but the shape of the eigenfunctions is sharply 
different). 

Namely, by Proposition 2.1, when j 7^ 1 and Xj-i,Xj are very near each 

other, the eigenvector V'^-'^ related to the eigenvalue \j : xj^i < Xj < xj has 

(7) (i) 

two main peaks of opposite sign given by V'j-i ^^'^ V'j j this is very different 
from the predictions of [24] (see their Figure 4). 

A. 3. Complex integral representation. Let L be a Markov generator on 
the state space S := {1, 2, . . . , A^}, reversible w.r.t. a positive measure /i. 
We can think of L as a linear operator on M^, symmetric w.r.t. the scalar 
product (-jOa*, where 

N 

{a,b)^ = ^fi{{)aibi. 
1=1 

In what follows we endow with the scalar product (and not 

with the standard Euclidean scalar product). Since L is symmetric, we 
can orthogonally decompose as = Wi © W2 © • • • © Wm such that 
L = J2k=i -^fc-FWfe) where Pw^. denotes the orthogonal projection on Wk and 
Aj / Xj if i / j. Given A G C \ {Ai, . . . , Am}, we write -R(A) for the resolvent 

m 1 

R{X):=iXI-Lr' = Y,j-y^Pw,. 
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Then, the residue theorem imphes the integral representation 
(A.ll) e"*^ = — / e-*^i?(A) dA, 

where 7 is a positive oriented loop containing in its interior Ai, A2, • • • , Am- 
Given a probability measure v on 5, we denote by the probability 
measure on the path space associated to the continuous-time random walk 
on S with generator L and initial distribution v. Fix j € S and let 
V G be such that Vi = 5ij. We write ^ for the Radon derivate, that is, 

Then the symmetry of L w.r.t. the scalar product (•, implies 



dfi 



N 

k=l 



By plugging (A.ll) in the r.h.s. of (A.12), we get the integral representation 

N 

(A.13) p,(Yit)=j) = -^ I e~'^{y2fiU)RjkW^\dX. 



In particular, given h function on S, 

Ep,(Ml^(t))) = EE^MjHj)^ J e~''Rjk{\)d\, 

J=l k=l ^ 

thus allowing to get an integral representation for Il{t, t^) ■= Pj/(no jumpin t«, + 
t]). If we set = Ti = and z^(i) = A^~^ (uniform initial probability), 
then 

(A.14) P„(y(*,=,-) = -L|^,-.A{^^|fl^,(„j,,. 

Let us consider now the special case given by Bouchaud's REM-like trap 
model where L := Ljv is defined in (2.3) and z> is the uniform distribution on 
5. Note that all the integral formulas obtained in Section 2 can be derived 
from the following one: 

(A.15) r„(rw=,) = ^/e-"(3^-i)^^A, 
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where (j){X) = J2k=i \-x ■ what fohows we prove that (A. 15) corresponds 
to (A. 14). 

We know aheady that det(AI — L) has distinct zeros given by the N 
distinct zeros of (j){X). In particular, it must be 

(A.16) det(AI-L) = l</>(A)n(A-x,) = ^A5] H (^"^^O- 

j k j:jy^k 

Given a matrix A, we write [A]ij for the determinant of the matrix obtained 
from A by erasing the ith. row and the jth column. Since 

^^'^^^^-^-^^ det(AI-L) 

and due to (A.16), in order to derive (A. 15) from (A. 14), we only have to 
show that 

(A.17) ^(_i).-+fc+i^[AI-L],,,= n (A-Xs). 

k s: Sy^j 

In order to prove the above identity, observe that [AI — ^]k,j is a polynomial 
of degree — 1 k = j, otherwise it has degree N — 2. The l.h.s. of (A.17) 
is a monomic polynomial of degree — 1. At this point, we only have to 
verify that Xs, s j, are zeros of the l.h.s. of (A.17). This is trivial if one 
observes that the l.h.s. of (A.17) is the determinant of the matrix obtained 
from AI — L by replacing the j'th column with the vector w with Wi = — 

Xj 

for i = 1, 2, . . . , A^. It is easy to verify that, if A = for some s / j, the jth 
row and the sth row in such a matrix are proportional, thus implying the 
thesis. 
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